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Abstract 

In this review a new method is presented for attaching parton shower algorithms to NLO 
partonic jet cross sections in electron-positron annihilation. Our method is based on the Catani- 
Seymour dipole subtraction method and also uses an adaptation of the matching scheme of Catani, 
Krauss, Kuhn, and Webber. 

1 Introduction 

One often uses perturbation theory to produce predictions for the results of particle physics ex- 
periments in which the strong interaction is involved. In order to get useable predictions one 
has to calculate at least at next-to-leading order to avoid large uncertainties those come from the 
unphysical scale dependences. 

Unfortunately, standard NLO programs have significant flaws. One flaw is that the final states 
consist just of a few partons, while in nature final states consist of many hadrons. A worse flaw is 
that the weights are often very large positive numbers or very large negative numbers. 

There is another class of calculational tools, the shower Monte Carlo event generators, such 
as HERWIG^ and PYTHIAEl. These have the significant advantage that the objects in the final 
state consist of hadrons. Furthermore, the weights are never large numbers. Finally, the programs 
have a lot of important structure of QCD built into them. For at least some cases like this, the 
shower Monte Carlo programs can provide a good approximation for the cross sections. The chief 
disadvantage of typical shower Monte Carlo event generators is that they are based on leading 
order perturbation theory for the basic hard process and thus reproduce only the first term in the 
perturbative expansion when applied to an infrared safe observable. 

It is possible to add the machinery of a shower Monte Carlo event generator to a next-to- 
leading order program in such a way that the complete program produces realistic final states 
made of hadrons, with weights are not unbounded in size. One example is the program of Frixione, 
Nason, and Webber^, which so far has been applied to cases with massless incoming partons but 
not to cases with massless final state partons at the Born level of calculation. The other example is 
that which concerns three-jet observables in electron-positron annihilation and thus addresses 
massless final state partons but not massless initial state partons. 

In this paper, we want, most of all, to have an algorithm that can be used by NLO practitioners 
in a reasonably straightforward manner. For this reason, we have based the algorithm on the 
dipole subtraction scheme of Catani and Seymour^. It is quite widely used for NLO calculations 
(for example in the programs NLOJET++® and MCFtM). 



We note that NLO calculations are generally limited to just one class of observables - for instance 
four-jet production but not at the same time three-jet production. We would like to overcome this 
limitation. For this reason we have adapted the fcr-jet matching scheme of Catani, Kuhn, Krauss, 
and WebbeJSlto the present circumstances. We also seek be as independent as possible of the choice 
of any particular shower Monte Carlo event generator. That is, we do not think that practitioners 
of NLO calculations should need to do separate calculations for each present and future shower 
Monte Carlo. 

The general idea of the algorithm that we present applies, we believe, to lcpton-lcpton collisions, 
lepton-hadron collisions, and hadron-hadron collisions. 

In the next section we have a very brief review of the algorithm. The precise definition and the 
details of the algorithm can be found in Ref.^l 

2 Structure of the algorithm 

The cross section computed with parton showers will consist of contributions from each available 
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For the contributions at NLO level, there are three terms, which correspond to Born, real emission, 
and virtual loop contributions with showers added ("+S"'). For the remaining terms there is only 
a Born contribution. We will arrange that (for a suitably behaved observable) 
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that the NLO expansion of the partial shower cross sections gives the correct NLO partonc cross 
section (<rJJ LO ) plus higher order and power corrections. 

2.1 Born term with showers 

Our discussion begins in this subsection with c^ +s . We define 
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The first line contains integrals over Born level parton momenta (dr({p} m )) and a corresponding 
sum over parton flavors. The second line contains sums over choices I of the parton that splits and 
a spectator parton k along with an integral over the splitting variable^ Yi — {yi, zi, 4>i, fi t \, /j 2} 
and a matrix element of certain operators Ei t k acting on the Born amplitude /} m )). The 

operators E^k together with the other factors in the formula describe the formation of showers 
from the Born level partons. 

The integration over the momenta {p} m is restricted by a factor 0(d m i < d m ({p, f} m )). Here 
d m ({p, f} m ) is defined by applying the kr jet finding algorithm to the m parton momenta. Given 
an n-parton final state, we apply the recursive "fcy" jet finding algorithm^" to the parton momenta 
{p, /}„, successively grouping the partons into jets. The algorthm gives a sequence of jet resolution 
parameters dj({p, /}„) at which two jets were joined, reducing J jets to J — 1 jets. Typically one 
has d n < d n -\ < ■ ■ ■ < d%. 

There is also a factor W m ({p, f} m )> which is the product of factors associated with the splitting 
history that matches the found jet structure, following the method of Ref.^. 



"The splitting variables are the yi virtuality like variable, zi momentum fraction variable, <f>i azimuthal angle and 
flavors of the emitted partons. The integral over these variables is 
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The function I({p, f} m ;l, k,Yi) is the interface function to secondary shower. Its important 
property is that the secondary shower provides only perturbative and power correction, thus we 
have 

I({pJ} m ;l,k,Y l )=F m+1 ({p} m+1 ) + 0(a s ) + 0(lGeV/V^) , (4) 

where the function -Fm+i({p}m+i) is the measurement function of an infrared safe observable. 

Now, we turn to the splitting function Ei ^ in Eq. @, which is an operator on the flavor and 
spin space of parton I in the vector \A4). This operator has the following form, 
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The parton splitting is organized according to an evolution parameter r, which is defined to be pro- 
portional to the transverse momentum square, Ri({p, f} m , y, z) = s; yz(\ — z) and s; is a virtuality 
scale appropriate to parton I. The simplest choice would be s; = s. . 

With the use of function 6 d({p, f} m , I, yi, zi), we limit the splitting in Ei to be unresolvable at a 
scale d that is approximately c?i n i x 2pi -phjsi. With this cut the vetoing procedure is implemented 
in the first step of the shower ensuring the cancellation of the c?i n i dependences al least at NLL 
level. 

In each E^f. operator, there is an operator on the parton color space, 3] • Tf], that is 

for the soft color connections. The splitting function Si acts in the spin space of the emitter and 
depends on the splitting parameters Y/ for parton I as well as on the momentum pi . This functions 
are proportional to the dipole splitting functions. 

The next factor, the Sudakov exponential, gives the probability that none of the partons has 
split at a higher evolution scale. The factor (S(yi>, Zv, /;/)) in the Sudakov exponent is the average 
over angle and flavors of S for parton I'. 

2.2 NLO corrections with shower 

We turn to the discussion of the NLO corrections. Let us start with the real contribution. Define 
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The first term is the m + 1 parton matrix element squared with the proper m-jet definition and 
the second term is the sum of the dipole contributions to eliminate the infrared singularities. The 
dipole function fy,fe are based on the m-parton tree level color connected matrix elements and the 
correct definition can be taken from the paper by Catani and Seymour^! 
The definition of the virtual correction is 



v+s 

yJ 71) 



^(W™) ( d «» < d rn({pJ}m)) I({p, f}m) 

x lv({ P J} m ) - W£\{p,f} m ) \M({ P J} m ) 



b lt is an approximation of the jet resolution variable, 
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where the function V({p, f} m ) represents sum of the 1-loop matrix element and integrated sub- 
traction term given in the second term of Eq. @ over the phase space of the unresolved particle. 
The function Wm* {{p, f} m ) is coefficient of the a s term in the expansion of Sudakov reweighting 
factor W m {{p, f} m )- 

2.3 Secondary shower 

Consider, the function I({p, f} m ; I, k, Yj) used for cr^ +s . In this term the E operator describe the 
emission of the hardest splitting, thus all the further splittings are constrained. 

All of the partons are allowed to split, and the one that does is parton / with aid of spectator 
k. The others did not split at an evolution variable above the value r. That is, parton /', with the 
aid of spectator k' , did not split with r' > r. Further evolution of these partons should then be 
restricted to the range r' < r. This constraint for the transverse momentums is 

A restriction like this can be imposed in the chosen shower Monte Carlo program by using a veto 
algorithm, as described for instance in RefM A sensible choice for the k' would be to let k! be one 
of the final state partons to which parton V is color connected (at leading order in 1/N C ). For the 
splitting of one of the daughters of parton I, one may simply impose \k±\ < \kj_\. 

Finally, in principle there should be a cut d(j>i,pj) < d ln i imposed on further splittings. However, 
for most events passed to the Monte Carlo, \kj_\ will be much smaller than d; n iS, so that this cut is 
not really needed. 

3 Conclusion 

We have proposed an algorithm for adding showers to next-to-leading order calculations for e + + 
er — > N jets. This algorithm is based on the dipole subtraction scheme^ that is widely used for 
next-to-leading order calculations. We also use the fey-jet matching scheme of Ref.^in order to 
incorporate the possibility of calculating infrared safe AT-jet cross sections for different values of N 
into the same computer program. 
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